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ELECTRON  SCATTERING  FROM  SILICON  30 


GEORGE  PAUL  BERNHARDT  IV 
ILt,  USAF 

1983 

61  pages 
S.M.  Physics 

MASS.  INSTITUTE  OF  TECHNOLOGY 
ABSTRACT 

An  electron  scattering  experiment  was  performed  at  the  MIT  Bates  Linear 
Accelerator  facility  on  a  target  of  n|gurally  occurring  silicon.  Several 
of  the  low-lying  inelastic  states  of  Si  were  examined  to  determine  the 
single-particle  transitions  which  contribute  to  the  states. 

The  levels  examined  where  the  second,  third,  eleventh,  twelth,  and  thir¬ 
teenth  levels.  Form  factors  for  each  of  these  states  were  determined 
experimentally  and  the  appropriate  harmonic  oscillator  polynomial  was  then 
extracted  from  the  form  factors.  Derived  polynomials,  using  a  harmonic  os¬ 
cillator  model,  for  single-particle  transitions  were  then  compared  to  the 
fitted  polynomials. 

The  second  and  twelth  levels,  J''  =  2^,  were  |round  to  be  mostly  a  2s-ld 
transition  while  the  third  level,  also  =  2_,  was  found  to  be  mostly  a 

Id-ld  transition.  The  eleventh  state,  =  3  ,  was  found  to  be^ approxi¬ 

mated  best  by  a  Id-lf  transition.  The^ghirteenth  level,  =  4  ,  led  to 
the  possibility  of  If  strength  in  the  si  ground  state  as  other  possible 
transitions  are  unlikely. 


CHAPTER  I 


INTRODUCTION 


In  1940,  the  first  experiments  on  electron  excitation  of  nuclei  to 
discrete  levels  were  done  by  Collins  and  Waldman(l. 1) .  The  first  theoreti¬ 
cal  discussions  of  this  inelastic  electron  scattering  were  due  to  Mamasach- 
lisor  in  1943(1.2)  and  Snedden  and  Touschek  in  1948(1.3).  Lyman  et  al.  in 
1951,  using  a  betatron,  clearly  saw  effects  due  to  deviations  from  point 
Coulomb  scattering,  thus  allowing  nuclear  sizes  to  be  measured  using  elec¬ 
tron  elastic  scattering(1.4) .  Since  that  time,  linear  accelerators  have 
become  the  standard  source  of  high-energy  electrons  for  such  work.  Hov/ever, 
due  to  the  limited  resolution,  an  experiment  of  the  type  done  in  this  thesis 
was  not  possible  before,  the  current  .generation  of  linear  accelerators  and 


'X.t-  \  0 

spectrometers.  ^  Jo 

The  machines  used  in  the  1950's  allowed  electron  scattering  with  a 


resolution  of  A  p/p  ^5x10  /and  with  currents  of  a  fe\<  tenths  of  a  micro¬ 
ampere.  Work  done  using  these  machines,  although  including  some  excellent 
work  such  as  that  done  by  Hofstadter (I . 5)  was  limited  to  a  few  nuclei  with 
well-isolated  levels  due  to  Inability  to  resolve  closely  spaced  levels.  By 

A- 3 

the  early  I970’s,  facilities  existed  with  energy  resolution  of  Ap/p  10 

with  useable  currents  of  l~2 /u{  A.  To  reach  a  resolution  on  the  order  of 
-4 

10  it  was  necessary  to  use  an  energy-loss  spectrometer. 


4 


In  previous  systems,  high  resolution  has  meant  throwing  away  most  of 

the  primary  beam  on  energy  defining  slits,  since  the  beam  incident  on  the 

target  was  analyzed  to  the  level  of  resolution  required  for  the  experiment. 

Less  than  one  percent  of  the  primary  beam  was  available  for  resolutions  on 
-4 

the  order  of  10  since  most  linear  accelerators  have  an  inherent  resolution 
-2 

on  the  order  of  10  .  The  most  important  improvement  in  electron  scatter¬ 
ing  of  the  last  decade  has  been  the  recognition  that  this  loss  of  intensity 
can  be  avoided  by  using  energy-loss  spectrometry.  The  two  important  reso¬ 
lutions  in  a  high  precision  measurement  of  an  electron  scattering  cross 
section  are  the  resolution  of  the  energy-loss  of  a  scattered  electron  re¬ 
flecting  excitation  of  nuclear  states  and  the  resolution  of  the  momentum 

transfer  to  the  nucleus.  Many  experiments,  including  this  one,  require  the 

-4 

nuclear  level  resolution  A  E/E  to  be  on  the  order  of  10  ,  but  the  defini¬ 

tion  of  raomentua  transfer  is  usually  adequate  at  ^q/q  on  the  order  of 

_3 

3x10  .  Energy  loss  spectrometry  exploits  this  relaxed  requirement  on 

momentum  transfer  resolution,  allowing  most  of  the  accelerator  beam  to  be 

used.  High  resolution  is  achieved  if  the  accelerator,  beam  transport  and 

spectrometer  form  an  achromatic  system  that  images  the  beam  on  the  focal 

plane.  Each  monochromatic  component  is  focused  on  the  target.  The  beam  is 

dispersed  spatially  at  the  target  to  match  the  spectrometer  dispersion(l .6) . 

-4  30 

Resolutions  of  the  order  of  10  are  required  to  perform  Si  scatter¬ 
ing  experiments  on  a  naturally  occurring  silicon  target.  Since  natural 
28 

silicon  is  92.2%  Si,  the  levels  from  this  isotope  will  be  dominant. 

30  28 

Therefore,  it  is  necessary  both  to  resolve  Si  levels  from  Si  levels 
and  ^^Sl  levels,  and  to  achieve  adequate  signal-to-noise  to  see  the  ^^Si 


levels  against  the  background  of  strong  radiative  tails  from  the  levels 


5 


This  thesis  presents  data  for  electron  scattering  from  several  of  the 
30 

low  lying  states  of  SI.  The  form  factors  for  these  levels  were  measured 
experimentally  and  compared  to  a  theoretical  form  factor  calculated  accord¬ 
ing  to  an  independent-particle  model  using  harmonic  oscillator  wave  func- 

28 

tlons.  Also,  the  form  factors  for  the  first  three  states  of  Si  were  com¬ 
pared  to  the  results  of  Whitner(1.7)  as  a  check  on  the  consistency  of  the 
data  analysis. 

Basically,  the  organization  of  this  thesis  is  given  below.  Chapter  II 
presents  the  theoretical  formulation  of  the  form  factors  using  harmonic 
oscillator  wave  functions  in  the  reduced  matrix.  Chapter  III  presents  the 
experimental  apparatus  and  procedure.  Chapter  IV  covers  the  data  analysis. 
Including  extraction  of  the  form  factors  from  the  spectra.  Chapter  V  com¬ 
pares  the  results  of  Chapters  II  and  IV  and  gives  conclusions. 
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CHAPTER  II 


ELECTRON  SCATTERING  THEORY 


INTRODUCTION 


The  Importance  of  electron  scattering  as  a  tool  for  studying  nuclear 
structure  cannot  be  understated.  The  electron  scattering  experiments  lead 
directly  to  an  understanding  of  the  electromagnetic  structure  of  nuclei. 
Three  considerations  show  this  Importance.  First,  the  electromagnetic 
interaction  between  the  electron  beam  and  the  target  is  well  understood, 
thus  allowing  separation  of  electromagnetic  effects  and  structural  effects. 
This  is  as  opposed  to  experiments  using  strongly  interacting  nuclear  pro¬ 
jectiles,  such  as  protons,  in  which  neither  the  interaction  nor  the  struc¬ 
ture  of  the  target  is  well  known.  Secondly,  within  the  restriction  that  the 
four  momentum  transfer  [q^  =(q,CJ)  ]  be  spacelike,  the  three  momentum  trans¬ 
ferred  to  the  nucleus  (q)  can  be  varied  independently  of  the  nuclear  excita¬ 
tion  energy,  allowing  one  to  map  the  transition  densities  to  discrete 


nuclear  states  in  momentum  space.  Thirdly,  in  many  experiments, 

X 


(1) 


is  small  compared  to  unity,  allowing  the  computation  of  the  transition  rate 
using  the  Born  approximation  for  exchange  of  one  virtual  photon.  In  such  a 
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case,  one  can  make  measurements  without  greatly  disturbing  the  structure  of 
the  target. 

This  section  presents  the  formalism  of  electron  scattering  largely 
based  on  the  development  by  Willey(2.1).  Emphasis  is  on  the  kinematics  and 
the  calculation  of  the  form  factor  using  a  harmonic  oscillator  model. 


KINEMATICS 


The  kinematic  parameters  are  shown  in  Figure  II. 1  for  an  electron- 
nucleus  scattering  event.  The  electron  is  extremely  relativistic  at  even 
the  lowest  Incident  energy  used  in  this  experiment  (148  MeV) .  Hence,  the 
electron  rest  mass  may  be  ignored  since  it  is  only  on  the  order  of  ,5  MeV 
Conservation  of  linear  momentum  yields  two  equations 

COS  9  i-  %  cos 

and 

O  =  SinO  +  %  Sin  (p  o 

Conservation  of  mass  and  energy,  in  a  relativistic  form,  yields 

Pi  Mt  •  Pf  *  ^  “ 


In  the  above  equations  Ex  is  the  excitation  energy  of  the  nucleus  and  Mj 
is  the  rest  mass  of  the  nucleus.  The  other  parameters  are  shown  in  the 
figure.  Using  these  equations,  one  may  derive  the  following  expression  for 


KINEMATIC  PARAMETERS 
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the  final  electron  momentuin: 


f>;  -  E«  -  ExVa  Mt 


(5) 


Cl  t 

^  (Hr 


J 


The  denominator  is  knovm  as  the  nuclear  recoil  factor  and  is  usually  denoted 
by  ^  .  If  one  solved  equations  2-4  for  q  Instead  of  for  P^,  the  results 


would  be 


where 


'  f/’:  Pf 


W  =  fc  -  Pi: 


This  leads  to  the  restriction  on  the  four  momentum  mentioned  earlier. 


(6) 


(7) 


~  >  O 


(8) 


In  analyzing  the  data,  the  kinematic  shift  caused  by  mass  differences 

between  Isotopes  must  be  taken  into  account.  The  excitation  energy  scales 

28 

for  the  spectra  are  determined  for  Si,  the  dominant  isotope.  The  peaks 
30 

due  to  levels  in  Si  will  appear  on  this  scale,  but  will  be  shifted.  The 


shift  arises  from  the  recoil  factor  because  of  the  mass  difference 
28  30  30 

between  Si  and  Si.  The  apparent  Si  levels  will  appear  to  be  lower 
than  their  true  energy.  If  Ex^  is  the  apparent  energy,  it  is  found  to  be 


(9) 


10 


where 


In  these  equations,  Ex2is  the  actual  excitation  energy.  The  convention  is 

28 

that  a  subscript  1  relates  the  factor  to  Si  and  a  subscript  2  is  the 


actual  value. 


SCATTERING  FORM  FACTOR 


Matrix  Formulation 

The  interaction  between  the  electron  as  it  scatters  and  the  nucleus 
leads  to  nuclear  transitions.  The  invariant  interaction  matrix  element  is 

<-f  J'^u) 


Here  A^(x) ,  the  electron  four-potential,  describes  the  electron  current  and 
is  composed  of  a  time-like  part,  the  scalar  notential  ^  and  a  space-like 
part,  the  vector  potential  A  (A^(x)  =  ( (x) ,A(x)) .  The  nuclear  four-current, 
J^(x) ,  is  composed  of  the  time-like  part,  the  nuclear  transition  charge 
density  (x) ,  and  the  space-like  part  composed  of  two  pieces:  (1)  j^^Cx), 
the  nuclear  current  transition  density,  and  {2)Vx  magnetiza- 

tion  transition  density  (J^(x)=  (x)^^;  cj^(x)j^  ) . 

For  a  gauge  invariant  interaction,  the  nuclear  current  must  be  con¬ 
served,  requiring  the  constraint 

^  ^  -  O  02) 
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The  Lorentz  gauge  places  the  constraint 


--  O 


(13) 


and  Maxwell's  equations  require 


D  -  -  e  Wf'  C  w.  <1 


e  _ 


(14) 


P 

where  D  is  the  d' Alembert ian  operator,  is  the  electron  four-current, 
is  the  appropriate  Dirac  matrix,  q^ is  the  four-momentum  transfer 
(q^  =  (q,  t*3)) ,  and  u^  and  u^  are  the  initial  and  final  electron  spinors.  The 
appropriate  Green's  function,  or  propagator,  for  this  equation  is,  by 


.1,  2, 


Bjorken  and  Drell(2.2),  -(  /q^)  yielding  the  solution 


f 
Uf' 


(15). 


The  transition  matrix  elements  can  now  be  written  as 


The  transition  probability  per  unit  time  is 


(16). 


t)  -  'Vs  / 


W( 


(17) 


where  is  the  spin  of  the  initial  state,  and  are  the  initial  and 
final  nuclear  spin  rates,  and  are  the  initial  and  final  electron  spin 

states,  and  (E^)  is  the  density  of  final  states  and  is  given  by 


where  is  the  final  electron  momentum,  dA^ is  the  element  of  solid  angle, 
and  is  the  recoil  factor  defined  in  equation  5.  The  cross  section  can 


now  be  written  as 


rfcr 

dJl 


J-  rj<f  1^-J_ 

3^^  \xXrj(^cJ  2.7^* t 


l<-flH'lL>l^ 


(18). 


The  cross  section  may  be  broken  up  into  components  consisting  of  the  coulomb 
multipole  term  (c) ,  the  transverse  electric  multipole  term  (E) ,  and  the 
transverse  magnetic  multipole  term  (M)  giving 

A*/  cA 


In  this  formulation,  transverse  and  logitudinal  components  are  defined  rela- 

u  . 

tive  to  the  direction  of  q.  In  the  above  equation, A  is  the  multipolarity 
of  the  transition  and  is  given  by 


I  If  -  J,  I  ^  ( 


Xf  *  Si 


(20). 


In  the  present  case,\=  since  J^*0  for  ^*^Si. 

The  transverse  multipoles  arise  from  a  spin  flip  of  the  incident  beam. 
Originally,  the  spins  of  the  electrons  are  aligned  in  either  the  direction 
of  the  electrons  (+)  or  opposite  the  direction  (-)  or  momentum.  After 
scattering,  the  electrons  must  conserve  momentum  and  spin.  The  final  state 
is  now  a  linear  combination  of  initial  +  and  -  states.  The  portion  which 
flipped  spin  from  either  plus  to  minus  or  minus  to  plus  accounts  for  the 
transverse  electric  and  magnetic  multipoles.  Thus,  the  longitudinal  coulomb 
multipole  results  from  that  portion  of  the  wave  function  which  did  not  flip 
spin. 

Just  as  important  as  the  origin  of  the  miltlpoles  is  the  parity  selec¬ 
tion  rules  governing  them.  The  coulomb  and  transverse  electric  multipoles 

A 

have  normal  parity  (-1)  ,  whereas  the  transverse  magnetic  multipole  has 
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abnomal  parity  (-1) 

Returning  to  the  formalism,  the  matrix  elements  from  equation  18  become 
<  -f  li>  =  ^ ^  ^  (^r)  \^(r)  ^ 

r aw  i)'/  c  -* 


.  uwi;r/  r  ^ 


(21b) 


KlU 

<flr^(MK-.i>)h>-<  t,; 


(21c) 


where  J  ^qr)  is  the  spherical  Bessel  function,  L  is  -ir  X"  V  ,  and  is  the 

spherical  harmonic.  These  equations  reduce  by  the  well-known  Wigner-Eckart 
theorem,  which  factors  out  the  dependence  on  the  magnetic  quantum  numbers  M. 


and  M^,  to  the  following: 


X  <  Tf  //  Hr>,{)// Ji> 


The  separated  cross  sections  may  now  be  written  as 

^  (]fCX  i  ^23a) 


14 


dA  ^  ' 'X^d^  Lott  *  ''''  a,  /  aCux^/;//J 
X  SCEX,i,^o-'>Jr) 


and 


^  \  X  fd<r  \  r  i  ^J.  CA»/)  j 

(JfXL  ~  '*1  '  dA  //vuTf  ^  /  ^  (.(’3lAv-/)//7^ 

X-  6  f  M  X ,  f  ,  r.  '?  ) 


where 


and 


X 


Ui^  c 

'i  cos  \ 

[3  £4 

s.„x  % 

/d£\ 

\d!A  l^atr 


and is  the  nuclear  recoil  term.  The  factor 


] 


t  / 

takes  out  the  lowest  order  of  q-dependence  of  B,  so  that 

f, « ^  6  <«.  > ,  f ,  r,-?  J, )  ^  (l^fr)  6.  ^  ■’ 


where  B^(M>^J^ “TJ^)  is  the  gamma-ray  reduced  transition  strength.  In 
notation  the  form  factors  may  be  written  as 

-  tJT  ,6(cX,^,T.^r. 

L( 


(23b) 


(23c) 


(24) 


(25) 


(26) 


(27) 


this 


(28a) 
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-- 

and 


( xjii  ] 

\  >  / 

V  6  ce 

%  ) 

X  13  ('m  A  , 


t  ,  T.- Jf) 


(28b) 


(28c). 


Individual  Particle  Model 

It  is  now  necessary  to  develop  equation  2A  by  using  an  appropriate 
model.  Before  applying  the  model,  it  is  instructive  to  briefly  examine  the 
general  case  as  presented  by  Donnelly  and  Walecka(2.3) .  In  second  quanti¬ 
zation,  the  multipole  operators  have  the  form 

&\  (a  ,  Al,  <«*•  /  Tp 

where  and  /3  form  a  complex  set  of  single  particle  quantum  numbers,  T  and 
refer  to  isospin  and  isospln  in  the  Z  direction, is  the  creation 
operator  for  the  state  and  is  the  destruction  operator  for  the 
state  .  Taking  the  matrix  elements  between  exact  nuclear  states  yields 

O  ,%)i  i>-  V 


where 


^1;  '■  \c:  c,ii> 


(31). 
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Equation  30  expresses  the  matrix  element  of  an  arbitrary  multipole  operator 
as  a  linear  combination  of  single-particle  matrix  elements  multiplied  by 
equation  31,  which  are  simple  numerical  coefficients.  The  basis  of  equa¬ 
tion  30  can  be  doubly  reduced  with  respect  to  both  angular  momentum  and 
isospin  yielding 


<f.  T^////W>  r  ^Vltl  7-,;^ 


(32). 


(f.) 


-  £  <a/// 


Once  the  set  of  numerical  coefficients  have  been  determined, 

the  exact  transition  matrix  element  of  any  multipole  operator  can  be  expressed 
as  a  linear  combination  of  single-particle  matrix  elements.  In  the  present 
case,  the  transition  was  assiimed  to  be  single-particle  only.  Therefore, 
was  assumed  to  be  a  delta  function. 

In  the  individual  particle  model,  the  nucleus  is  considered  to  be  a 
collection  of  nucleons.  Therefore,  the  nuclear  charge,  current,  and  mag¬ 
netization  density  operators  are,  in  Willey's  notation, 

■  1^1  £j  r- 

‘'’’r  '■  'ii  •?,  hy... 


(33a) 

(33b) 

(33c) 


where 

(S.  =  -  ,  V. 


(34). 


The  symmetrized  form  for  equation  33b  is  necessary  for  the  nuclear  current 
in  order  to  satisfy  hermiticity  and  the  continuity  equation 


i  < 


(35). 
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in  isobaric  spin  notation 

i  0  <  ’  4  )  f., 

y.  z  \  (i  *  ■*  i  *'•' 


(36a) 


(36b) 


2.79  and  =  -1.91.  Using 


where  "C^  -  2t^,  =  It  ■  0. 

these  results t  the  inultipole  operators  become 

H^,>)  Y,^  (r,)  +  J  ( 

+  3>,  ^  ^  »,  if  Wb) 

.  V  /  ^  \  .  M  .■  \  .  _ 


X  n ) . /.  ^  )./» ^  ^ ^  ^ 

*  >'i  t- i  »  (  TO  )'''■  j,..  >  V,,„^  Cr,  )  •?; 


(37c) 


In  the  above  equations,  j(4rj)  is  the  spherical  Bessel  function,  Yj^^^(tj) 
is  the  vector  spherical  harmonics,  Y^^(r^)  is  the  spherical  harmonics, 
and  i;ls  -Ir^X^-. 

The  last  step  in  developing  the  form  factors  is  to  use  a  model  for 
the  initial  and  final  states  of  equation  24.  In  this  case,  a  harmonic 
oscillator  model  was  used.  The  reduced  matrices  have  been  calculated  by 
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Donnelly  and  Haxton  using  this  model.  The  general  result  Is 

i3  ,  r,  ^ ->  y  V  p  ’'‘r> 

where  ^  can  be  related  to  the  gamma  ray  transition  strength,  y  Is  mK 
b  Is  the  harmonic  oscillator  length  parameter,  and  p(y)  Is  a  polynomial 
of  y  whose  coefficients  are  given  In  the  tables  by  Donnelly  and 
Haxton(2 .4) . 

In  this  work,  the  transverse  electric  and  magnetic  multipoles  are 


assumed  to  be  very  small  for  the  data  taken.  Hence,  only  the  longitudinal 
coulomb  form  factor  Is  calculated  In  Chapter  V. 
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CHAPTER  III 


EXPERIMENTAL  APPARATUS 


INTRODUCTION 


The  data  presented  in  this  thesis  are  the  results  of  an  electron¬ 
scattering  experiment  on  natural  silicon  performed  at  the  MIT  Bates 
Linear  Accel-'rator .  The  discussion  in  this  chapter  is  limited  to  a  des¬ 
cription  of  those  unique  features  of  the  facility  relevant  to  this  exper¬ 
iment,  as  the  accelerator  has  been  described  in  detail  in  the  literature 
(3. 1-3. 2), 

The  facility  provides  an  average  electron  beam  current  of  2-40 A 

without  sacrificing  resolution  (^2  x  10  ^  through  its  unique  energy- 

loss  spectrometer.  The  high  beam  current  was  necessary  to  provide  a 

reasonable  measurement  of  the  small  experimental  cross  sections  due  to 

30 

the  low  percentage  of  Si  in  the  target.  In  natural  silicon,  the  percent¬ 
age  of  ^^Si  is  only  3.1%.  The  resolution  must  also  be  high  enough  to  pro¬ 
vide  a  reasonable  slgnal-to-noise  ratio  in  the  measured  spectrum.  Without 
the  high  beam  current  and  good  resolution,  the  experiment  would  have  been 
impossible  with  natural  silicon. 
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SPECTROMETER  AND  DATA  ACQUISITION 

In  the  Bates  energy-loss  system,  the  Incident  beam  dispersal  on  the 
target  allows  the  position  on  the  target  to  be  correlated  with  the  incident 

4 

energy  to  an  accuracy  of  better  than  one  part  in  10  .  A  quadrupole  sing¬ 
let  controls  the  dispersion  at  the  target  in  order  to  match  the  dispersion 

of  the  beam  transport  system  to  that  of  the  spectrometer.  In  order  to 
-4 

achieve  10  resolution,  the  scattered  electrons  are  bent  through  ninety 
degrees  by  the  split-pole  spectrometer  in  proportion  to  their  momentum 
and,  therefore,  their  position.  This  differential  bending  brings  electrons 
that  underwent  the  same  energy  loss  in  the  target  to  a  common  focus  in 
the  spectrometer  focal  plane(3.3). 

In  the  detector  system,  an  electron  passes  first  through  two  multi¬ 
wire  proportional  counters,  then  into  two  lucite  Cerenkov  counters.  The 
planar  surfaces  of  all  four  detectors  are  at  an  angle  of  forty-five 
degrees  to  the  central  ray  of  the  spectrometer. 

The  first  proportional  counter  has  sense  wires  tapped  by  three  differ¬ 
ent  delay  lines  read  by  tlme-to-digital  signal  convertors.  The  drift  time 
of  the  ions  to  these  wires  can  be  used  to  measure  both  the  track  position 
and  the  angle  of  entry  into  the  chamber.  Drift  times  are  measured  perpen¬ 
dicular  to  the  plane  of  the  sense  wires  and,  therefore,  this  counter  is 
known  as  the  vertical  drift  chamber.  The  sense  wires  give  information 
about  particle  trajectories  in  the  direction  of  particle  momentum. 

The  second  proportional  counter  is  known  as  the  transverse  array  and 
has  sense  wires  read  on  only  one  delay  line  in  the  direction  transverse 
to  the  momentum  of  the  particle.  Since  the  optics  of  the  spectrometer  are 
point-to-parallel  in  this  plane,  this  chamber  provides  a  measurement  of 
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the  scattering  angle  of  the  particle  within  the  acceptance  of  the  spec¬ 
trometer  (3.  4)  • 

The  coincidence  between  the  two  Cerenkov  counters  provides  a  start 
signal  for  the  delay  time  measurements  of  the  tlme-to-dlgltal  signal  con¬ 
vertors.  A  coincidence  between  two  Cerenkov  counters  Is  used  due  to  the 
need  for  efficiency.  The  discriminator  on  each  Cerenkov  must  be  set  very 
low  for  good  efficiency,  which  would  result  in  a  high  noise  background  if 
only  one  were  used.  The  Cerenkov  threshold  also  discriminates  against 
heavy  charged  particles  whereas  a  plastic  scintillator  would  not. 

Storage  of  the  event  parameters  (momentum  position,  vertical  angle, 
and  horizontal  angle)  and  event  classification  are  done  by  the  data 
acquisition  programs  in  the  on-line  PDP-11  computer.  Communication  be¬ 
tween  the  PDP-11  and  the  detectors  is  done  through  the  Caraac  Dataway 
interfaced  by  a  microprogrammed  branch  driver. 

The  detection  system  requires  300  nano  seconds  to  process  the  infor¬ 
mation  from  a  given  event.  Should  another  event  occur  during  this  time 
span,  both  the  start  and  the  new  event  are  disqualified  by  the  system  and 
stored  in  a  special  register.  The  recorded  signals  are  tested  by  known 
timing  criteria  to  determine  whether  the  signals  are  consistent  with  a 
"good"  event.  Inconsistent  signals  are  also  disqualified  and  stored  in 
registers  according  to  the  criterion  covering  the  disqualification. 

Typical  resolution  values  for  this  experiment  ranged  from  1.6  x 

10  ^  to  2.0  X  10  ^  inside  the  range  necessary  for  ^®Si.  Average  beam 

P 

currents  on  the  target  during  the  experiment  ranged  from  4  to  40  a^.A, 
with  20-30  most  typical. 
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TARGETS 


The  targets  used  in  this  experiment  were  natural  silicon  crystals 
provided  by  Or.  Robert  Wolf son  and  Mr.  Paul  Grubenskas  of  Varian  Vacuum 
Division.  The  targets  were  etched  to  their  thickness  used  in  the  experi¬ 
ment  from  10  mil  thick  blanks  using  a  blend  of  hydrofluoric,  nitric,  and 
acetic  acids.  No  contamination  as  a  result  of  this  process  was  observed. 

A  carbon  elastic  peak  was  observed  in  three  of  the  data  runs  used  for 
this  analysis.  However,  this  contamination  was  most  likely  an  impurity 
in  the  vacuum  system,  and  not  a  result  of  the  etching  process.  The  car¬ 
bon  peaks  did  not  interfere  with  any  of  the  data  analysis.  No  other  con¬ 
tamination  was  observed  at  the  level  of  this  experiment. 

An  analysis  of  target  thickness  was  done  by  Whitner(3.5) .  Summari¬ 
zing  her  work,  she  found  the  thickness  of  the  standard  target  by  using  the 
focal  plane  calibration  program,  FPCAL  1,  described  later.  The  silicon 
target  was  used  with  other  targets  of  known  thickness  to  calculate  a  fit 
to  the  focal  plane  and  accelerator  energy  parameters  of  the  run.  A  map 
of  the  fit  chi-square  as  a  function  of  target  thickness  was  produced  by 
varying  the  silicon  target  thickness  in  the  calibration  program.  The 
target  thickness  corresponding  to  the  minimum  of  the  chi-square  was  chosen 
as  correct. 

Target  thicknesses  were  compared  by  measuring  the  form  factor  of  the 

28 

Si  2+  level  under  the  same  conditions  with  the  four  individual  targets. 
The  ratio  of  these  form  factors  to  the  standard  target  form  factor  was 
then  used  to  find  the  other  target  thicknesses.  A  summary  of  the  thick¬ 
nesses  according  to  the  measurements  of  Whltner  is  shown  in  Table  Ill.l. 
Target  uniformity  was  measured  by  attenuation  of  225  MeV  electrons  from 
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Pm  in  the  macerlal(3.6) .  The  average  thickness  in  the  beam  spot  area 
was  found  to  vary  -3%  In  the  center  of  the  target. 


Table  (3.1) 
Target  Thicknesses 


Target  Thickness 
(mg/cm^) 


26.50 


25.37 

28.83 


31.83 


Uncertainty 


24 


CHAPTER  IV 


DATA  ANALYSIS 


INTRODUCTION 


A  raw  experimental  datum  is  simply  the  number  of  events  per  detector 
channel.  To  relate  these  events  to  the  form  factors  of  Chapter  II,  the 
relationship  between  the  experimental  parameters  and  those  scattering 
properties  of  interest  must  be  used  to  reduce  the  data  to  useable  form. 

The  reduction  is  done  by  creating  a  spectrum  of  true  events  as  a  function 
of  nuclear  excitation  energy,  by  finding  the  area  of  the  peaks  in  that 
spectrum,  and  by  then  calculating  the  form  factors  of  those  peaks  from 
that  area  using  the  known  values  of  target  thickness  and  integrated  change. 


ENERGY  SPECTRUM 


In  order  to  determine  the  nuclear  excitation  energy  corresponding  to 
a  given  detector  channel,  it  is  necessary  to  use  the  design  parametrization 
of  the  spectrometer  and  the  focal  plane.  In  this  experiment,  two  programs 
were  used  to  achieve  the  spectrum,  FPCAL  1,  which  determines  the  focal 
plane  parameters  and  accelerator  energy,  and  SORT,  which  combines  spectra 
at  the  same  magnet  settings. 


FPCAL  1,  although  used  for  some  time  at  MIT,  has  not  been  documented 
to  the  extent  given  here.  The  model  used  in  the  analysis  is  Model  1,  the 


model  used  for  the  design  of  the  spectrometer.  In  general,  the  program 
performs  a  least  squares  fit  to  determine  three  focal  plane  parameters 
and  the  bombarding  energy. 

The  kinematics  are  those  of  Chapter  II  with  one  addition.  The  rela¬ 
tionship  between  the  scattered  electron  momentum  and  incident  momentum  is 
given  by 


f-  E,  -  f//aM 


i  + 


4 


C39) 


where  ^  is  the  most  probable  energy  loss  in  the  target  and  all  other 
parameters  are  the  same  as  in  Chapter  II.  The  quantity  A  is  given  by 
Landau (4.1)  as 

A  -  S  O’  +  ^‘'3  J  ^  ^  /^eff  )  ]  ( 

where 

and 

2,ff  - 


(41) 


(42) 


The  quantities  in  the  above  equations  are  defined  as  f^^  being  the  fraction¬ 
al  abundance  of  the  i'th  species  in  the  target,  being  the  atomic  num¬ 
ber  of  the  i'th  species  in  the  target,  being  the  atomic  weight  of  the 

i'th  species  in  the  target,  being  the  angle  between  the  beam  and  the 

2 

normal  to  the  target,  and  a  is  the  target  thickness  in  mg/ cm  .  The 
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unknovm  quantity  which  will  be  fitted  from  equation  39  is  the  incident 
momentum,  P^. 

The  displacement  from  the  central  ray  along  the  focal  plane  is 
assumed  to  be 


(43). 


Defining  £,  to  be  X ‘f/^ ),  yields 


(44) 


where  x  is  the  distance  from  the  central  ray,  (x/^)^  is  the  first  order 
dispersion,  /  is  the  second  order  dispersion,  is  defined  above,  and 
is  the  electron  momentum  corresponding  to  the  central  ray.  The  form  of 


P^  is  given  by 


(43) 


where  is  a  linear  constant  with  a  derign  value  of  66.90  MeV/KG,  and  B 

is  the  magnetic  field  strength  in  KG.  The  form  of  P^  may  be  revised  to 
include  possible  non-linearities  if  necessary.  The  unknowns  from  equation 
44  which  must  be  fitted  are  the  two  dispersion  terms.  Equation  4^  may 
be  solved  for  P^  to  yield 


f’f 


<f6.y 


(46). 


This  equation  may  be  expanded  and  approximated  by 


(47). 


For  the  Bates  spectrometer ,  6, 'S  .4,(>r/i)  S  372  inches  and  X  F  12  in. 
These  terms  are  good  initial  guesses  for  the  least  squares  fit.  Using 
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r  - 


4 

these  approximations,  equation  47  should  be  valid  to  about  1  x  15  over 
the  entire  focal  plane. 

The  data  acquisition  code  used  in  this  experiment  creates,  by  soft¬ 
ware,  3072  channels  of  equal  space  across  the  focal  plane  into  which  the 
data  are  collected.  The  distance  X  can  be  obtained  from  the  focal  plane 
geometry  by  the  equation 


y  ^  c/f  54  - 


7)W  -  ( 


t  ooc 


) 


(48) 


where  1536  is  the  central  channel,  J  is  the  channel  in  question,  W  is  the 
width  of  the  channels,  presently  1/128  in.,  F  is  the  focal  plane  reference 
reading,  presently  1500,  and  Z  is  the  readout  value  when  the  focal  plane 
array  is  positioned  so  that  the  central  ray  passes  through  the  central 
channel  of  the  focal  plane  array.  In  equation  48,  the  quantity  Z  is  the 
quantity  which  must  be  fit. 

To  calibrate  the  energy  spectrum,  one  must  use  known  peaks  in  the 
spectrum  to  determine  the  focal  plane  parameters.  Combining  equations  44 
and  48  and  solving  for  the  channel  number  in  which  the  i'th  peak  falls, 
yields 


The  energy  of  the  peak  must  be  known  and  associated  with  the  i’th  channel. 
The  program  now  has  a  series  of  known  J^'s  associated  with  a  number  of 
known  peaks.  A  least  squares  fit  is  now  performed  on  equation  49,  varying 
the  parameters  P^,  Z,  X/f ^  and  6,  .  The  quantity  is  related  to  the 
quantity  P^  in  equation  49  by  equation  39. 


M 
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The  program  FPCAL  1  can  also  cake  into  account  the  change  in  momentum 
from  one  end  of  the  focal  plane  to  the  other.  In  Model  2,  equation  43 
becomes 


fit 


(50). 


Solving  this  equation  for  P, 

P  - 

f  aA 


yields 

-  ? 


(51). 


In  equation  51,  the  first  terra  explicitly  cakes  into  account  the  change 
in  momentum  from  one  end  of  the  focal  plane  Co  Che  ocher  and  Che  second 
term  corrects  for  spatial  variation  in  this  change  in  width.  The  procedure 
from  this  point  on  is  identical  to  that  of  Model  1. 

The  program  SORT  uses  the  focal  plane  parameters  to  produce  a  spec¬ 
trum  of  events  as  a  function  of  excitation  energy.  The  spectra  can  be 
reorganized  into  channels  of  uniform  excitation  energy  size.  Since  each 
detector  channel  has  a  spread  Ap/p  •  2  x  10  the  energy  size,  iip,  depends 
on  the  value  of  p.  A  rehistogramming  subroutine  in  the  program  SORT  will 
reorganize  the  data  into  bins  of  uniform  excitation  energy  width. 

The  program  SORT  will  also  combine  data  files  taken  at  the  same  inci¬ 
dent  energy,  scattering  angle,  and  magnet  setting,  but  at  different  focal 
plane  positions.  The  distance  X  along  the  focal  array  is  given  in  terms 
of  a  known  reference  position.  Therefore,  data  from  another  focal  plane 
position  can  be  added  by  specifying  the  new  reference  position. 

The  program  SORT  also  uses  information  about  the  total  charge  accumu¬ 


lated  in  each  data  run  by  the  beam  charge  integrators.  The  integrator 
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calibration  constant  Is  used  with  this  Information  to  determine  the  events 
as  counts  per  microcoulomb  of  Incident  charge. 

The  final  sorted  data  file  contains  the  events  In  counts  per  micro- 
coulomb  as  a  function  of  nuclear  excitation  energy.  The  energy  channels 
have  been  reorganized  Into  a  constant  width  across  the  focal  plane.  A 
typical  sorted  spectrum  Is  shown  In  Figure  IV. 1. 


CROSS  SECTION 


The  program  FITAE  was  used  to  extract  the  cross  sections  from  the 
sorted  energy  spectra.  The  code  uses  a  least  squares  fitting  routine,  by 
J.  Bergstrom,  to  approximate  a  line  shape  to  the  peak  in  question.  By 
integrating  over  some  energy  range,  the  peak  area  can  be  determined  and, 
thus,  the  differential  cross  section  may  be  determined  as  It  is  propor¬ 
tional  to  the  peak  area. 

The  shape  of  the  peak  can  be  fit  with  a  deformed  Gaussian  function  of 
the  form 

)  (52) 

where  HT  is  the  height  of  the  peak,  R  and  P  are  parameters  to  be  fit,  and 

~  ^  )  /  u/  >  *■( -h  l\ 

where  E  is  the  excitation  energy  (MeV) ,  E  peak  is  the  excitation  energy  of 
the  peak,  and  width  is  the  full  width  of  half  maximum.  The  deformed  Gauss¬ 
ian  is  joined  to  an  inverse  polynomial  tall  of  the  form 


Anti  IN 


i 


L_ 
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I  A  J  -  A  /x  I  t  6/yi^  ^  c/xr 

where  XI  >  X-S.  These  radiative  calls,  which  appear  on  the  lower  scattered 
electron  energy  side,  are  mainly  due  to  electrons  that  have  emitted  a 
photon  of  sufficient  energy  to  reduce  their  final  energy  below  the  energy 
cutoff.  The  peaks.  In  the  program,  are  defined  to  this  cutoff  energy, 
given  by  ^  E,  which  Includes  parts  of  the  peak  and  radiative  tail.  Radia¬ 
tion  calls  of  preceding  peaks  form  a  background  to  each  succeeding  inelastic 
peak  in  the  spectrum.  The  parameters  of  the  tail  function,  A  and  B,  are 
evaluated  internally  by  the  program  from  continuity  conditions  at  the  join¬ 
ing  point.  The  radiative  tail  is  joined  to  the  peak  function  at  the  point 
YM,  where  the  total  function  and  its  first  derivatives  must  be  continuous, 
and  YM  is  given  by 

YAl  '  e  (55) 

where  E'  is  the  excitation  energy  at  the  joining  point. 

The  parameters  HT,  E  peak,  width,  R,  P,  S,  C,  and  YM  can  be  varied  to 
minimize  the  chi-square  of  the  fit.  In  practice,  R,  P,  S,  C,  and  YM  are 
varied  for  a  reference  peak  only.  These  same  parameters  for  all  other 
peaks  are  then  fixed  equal  to  those  of  the  reference  peak.  This  fact  en¬ 
sures  that  peak  shapes  do  not  vary  widely  from  peak  to  peak.  Initial 
guesses  for  HT  and  width  were  obtained  by  examination  of  the  sorted  spec¬ 
trum.  The  initial  guess  for  E  peak  was  the  known  excitation  energy  for 
that  peak.  The  initial  guesses  for  the  other  five  parameters  were  the  de¬ 
fault  values  which  have  been  sec  by  experience. 
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Background  events.  Including  the  radiative  tails  mentioned  above, 
must  also  be  fit  to  Insure  proper  subtraction.  The  model  used  in  this 
experiment  is  given  by 

-  i\cn  T  ^  2^  (56) 

^  =  £  -  Ai*') 

where  the  four  A  parameters  can  also  be  varied  to  minimize  the  chi-square 
fit.  An  inverse  polynomial  background  can  also  be  chosen  but  was  not 
used  in  this  work. 

The  errors  associated  with  the  fitting  technique  were  also  calculated 
by  the  program  according  to  the  algorithm  established  by  Bevington(4 . 2) 
using  statistical  weighting  factors  and  the  inverse  derivative  matrix. 

The  peak  function  and  radiative  tail  are  integrated  to  the  cutoff 
point  as  mentioned  above.  At  this  point,  theoretical  radiation  corrections 
can  be  applied  to  obtain  the  total  peak  area(4. 3-4 . 4) .  These  corrections 
are  documented  fully  by  Whitner(4.5)  and  will  be  summarized  here.  The 
total  correction  to  the  experimental  cross  section  can  be  represented  by 
a  multiplicative  factor  C_  by 

ir  \  ^ (57) 

d-O. 

where  is  given  by  Fj^  F2  F^-  F^  is  the  Schwinger  factor  and  is  the  most 
important  correction.  The  term  arises  from  the  coherent  emission  and 
reabsorption  of  virtual  photons  by  the  electron,  and  the  coherent  emission 
of  real  photons  during  the  scattering  event.  As  stated  above,  the  majority 
of  the  radiative  tail  is  from  electrons  that  have  emitted  a  photon  of 
sufficient  energy  to  reduce  their  final  energy  below  ^E.  The  term  F^  is 
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the  Brensstrahlung  factor  and  arises  because  the  electron  passes  through 
a  target  of  finite  thickness.  Hence,  the  electron  can  be  deflected 
slightly  by  other  nuclei  before  and  after  the  scattering  event.  Thus, 

F2  Is  a  "thick  target"  or  Incoherent  bremsstrahlung  correction.  Finally, 
F^  Is  the  Landau  straggling  factor.  This  factor  corrects  for  energy  loss 
straggling  which  occurs  when  the  electron  Ionizes  an  atom  In  the  target. 
The  total  peak  area,  after  the  corrections  have  been  applied,  should  be 
Independent  of  the  cutoff  energy  ^  E. 

To  ensure  reasonable  values  for  peak  areas,  each  fit  was  evaluated 
on  whether  the  reduced  chi-square  was  small,  generally  less  than  3,  and 
If  the  corrected  peak  areas  were  independent  of  /lE  to  better  than  3%. 

This  last  step  was  done  by  using  five  different  cutoff  energies.  Finally, 
when  very  small  counting  statistic  uncertainties  meant  a  large  reduced 
chi-square,  the  human  eye  was  the  final  judge  of  the  acceptability  of  the 
fit. 

After  determining  the  peak  area,  the  differential  cross  section  was 
calculated  by  the  equation 

££  .  'jj  /  -1-  \  (38) 

dSL  "  ^  A  / 

Yq 

where  ^  is  the  peak  area  in  counts/microcoulomb,  e  is  the  change  of  one 
electron,  ^  -0.  is  the  spectrometer  solid  angle  acceptance,  and  is  the 
number  of  nuclei  in  the  target  per  unit  area. 

The  spectrometer  solid  angle  is  by  horizontal  and  vertical  slits  set 
by  the  experimenter.  Most  of  the  work  in  this  case  was  done  at  a  Viorizon- 
tal  setting  of  two  inches  and  a  vertical  setting  of  ten  inches,  giving  a 
solid  angle  of  3.325  millisteradlans . 
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The  number  of  target  nuclei  can  be  calculated  from 


A',  - 

w 


(59) 


where  Is  Avogadro's  number,  tg££  is  the  effective  target  thickness  in 
2 

mg/cm  ,  N  is  the  number  of  atoms  of  species  per  molecule,  and  U  is  the 
molecular  weight  of  species  in  mg/moles.  The  effective  target  thickness 
depends  on  the  orientation  of  the  target  with  respect  to  the  beam.  The 
data  were  taken  at  scattering  angles  of  45  and  90  degrees.  This  fact 
allows  the  targets  to  be  oriented  in  a  transmission  mode  in  the  target 
chamber  (Figure  IV. 2).  In  this  mode,  the  effective  thickness  of  the  tar¬ 
get  to  an  electron  is  independent  of  where  the  scattering  occurs  in  the 
target.  An  electron  scattering  anywhere  in  the  target  travels  through  the 
same  amount  of  material.  The  effective  thickness  is  given  by 


t 


^/c 


e>/3. 


(60) 


where  t  and  ^  are  defined  in  Figure  IV. X 


FORM  FACTOR 


The  form  factor  may  then  be  calculated  for  the  differential  cross 
section.  The  multipole  of  interest  was  the  coulomb  raultipole.  Since  only 
forward  angles  were  used  to  take  the  data,  the  coulomb  term  will  be  much 
greater  than  the  transverse  terms.  Combining  equations  23a,  24,  and  26 
from  Chapter  II  yields  a  differential  cross  section  of  the  form 


is" 

d£L  " 


-i-  (  ) 


IFI^ 


(61) 
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Therefore,  the  form  factor  squared  may  be  determined  by  multiplying  the 
differential  cross  section  by  the  nuclear  recoil  factor  and  by  Z  ,  then 
dividing  by  the  Mott  cross  section.  The  program  FITAB  calculated  both  the 
differential  cross  sections  and  the  form  factors. 

The  polynomial  of  equation  38  was  then  extracted  by  using  the  program 
SMFIT.  The  program  uses  a  least  squares  fitting  routine  to  determine  the 
coefficients  of  the  polynomial.  The  parameters  of  the  program  are  the 
harmonic  oscillator  length  parameter  of  Chapter  II,  the  single  particle 
normalizing  factor,  and  the  coefficients  of  the  polynomial. 

The  experimental  form  factors  for  a  given  transition  were  divided 

by  the  terms  in  front  of  the  polynomial  in  equation  38.  The  remaining 

function  was  then  fitted  by  a  least  squares  fit  to  determine  the  form  of 

the  polynomial.  The  results  of  the  fit  could  then  be  compared  directly  to 

the  tables  of  Donnelly  and  Haxton  to  determine  if  single  particle  transi- 

30 

tions  are  realistic  approximations  to  those  that  occur  in  Si. 

The  results  of  the  data  analysis  are  given  in  Chapter  V  and  discussed 
in  detail  there. 

SYSTEMATIC  ERRORS 

In  addition  to  the  statistical  errors  from  the  program  FITAB,  there 
are  several  sources  of  systematic  errors  listed  in  Table  4.1  and  discussed 
below. 

Errors  from  scattering  angle  measurement,  change  integration,  detec¬ 
tor  inefficiency,  and  solid  angle  measurement  are  errors  t  ultlng  from 
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the  limitations  in  the  actual  physical  measurement  of  equipment  parameters, 
errors  in  energy  calibration,  count-rate  corrections,  and  target  parameters 
are,  in  principle,  difficult  to  generalize  and  must  be  determined  separ¬ 
ately  for  each  data  run.  However,  in  practice  the  main  limitation  on  ex¬ 
perimental  accuracy  is  the  precise  knowledge  of  the  target  thicknesses 
and  uniformities,  while  other  corrections  vary  little  from  run  to  run. 

Table  4.2  shows  the  different  energies,  angles,  and  targets  used  to  provide 
the  data  in  this  experiment. 

The  energy  calibration  uncertainties  depend  on  the  individual  data 
run.  The  error  in  machine  energy  measurement  was  less  than  .2%  for  most 
of  the  runs. 

The  uncertainties  in  the  target  parameters  were  determined  by  Whitner 
and  form  the  largest  contribution  to  the  systematic  errors.  The  target 
thickness  was  known  to  within  5%.  The  error  associated  with  assuming 
natural  abundances  in  the  target  was  probably  less  than  1%. 

It  is  difficult  to  combine  systematic  errors  from  different  sources 
because  they  do  not  obey  statistical  distribution  laws.  However,  the 
best  upper  limit  estimate  of  the  systematic  errors  for  each  data  point  is 
7-10%. 
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Table  4.1 

Main  Sources  of  Systematic  Errors 


Machine  Energy  Calibration 

-  .3  MeV 

Scattering  Angle 

-  .05  degrees 

Charge  Integration 

1.0% 

Solid  Angle 

.5% 

Overall  Detector  Inefficiency 

2.0% 

Dead-time  Corrections 

.20  (f) 

Target  Thickness 

5-10% 

Target  Composition 

1.0% 

Target  Uniformity 

3.0% 

i 
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Table  4.2 


Incident  Energy 

Angle  (degrees) 

Target  Thickness 

148.678 

45 

26.5 

174.503 

45 

26.5 

199.141 

45 

26.5 

224.552 

45 

26.5 

166.360 

90 

26,5 

170.814 

90 

26.5 

199.994 

90 

31.83 

228. 780 

90 

26.5 

251.036 

90 

28.83 

265.451 

90 

28.83 

279.569 

90 

26.5 

292.968 


90 


26.5 
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CHAPTER  V 


RESULTS  AND  CONCLUSIONS 

This  chapter  will  examine  each  state,  comparing  the  extracted  polynomial 

with  the  polynomials  calculated  by  Donnelly  and  Haxton(5.1)  using  a  harmonic 

oscillator  model.  Also,  the  form  factors  calculated  for  the  first  three 
28 

states  of  Si  will  be  compared  to  those  of  Whitner(5.2)  to  ensure  there  are 
no  anomalies  in  the  data. 

The  form  factors  extracted  from  the  data  by  the  program  FITAB  are  shown 
in  Tables  V.l  through  V.6  and  in  Figures  V. 1  through  V.6.  Each  figure  is  di¬ 
vided  into  two  parts.  The  upper  portion  is  the  form  factor  plotted  versus 
q  effective  and  the  lower  portion  is  the  extracted  polynomial  of  Equation  38, 
which  is  given  as  p(y).  It  is  this  polynomial  which  may  be  compared  to 
those  derived  by  Donnelly  and  Haxton.  The  form  factor  plot  contains  three 
sets  of  data,  the  experimental  form  factors,  the  form  factors  calculated  from 
the  fitted  polynomial,  and  the  form  factors  determined  in  a  previous  electron 
scattering  experiment  by  Brain,  et  al(3.3). 

The  levels  shown  are  levels  two,  three,  eight,  eleven,  twelve,  and  thir¬ 
teen.  The  eighth  state  is  an  unnatural  parity  state  and  is  transverse  mag- 
netic.  Hence,  the  polynomial,  p,  cannot  be  extracted  in  the  same  manner  and 

is  not  included.  Other  levels  are  not  included  for  two  reasons.  First,  for 
30  28 

some  Si  levels,  the  Si  peaks  lie  too  close  in  the  spectrum  to  allow  mean- 

30 

Ingful  separation.  Secondly,  some  Si  peaks  are  so  weakly  excited  that  they 
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TABLE  V.  1 


Level  Two,  =  2  Excitation  Energy  2.234  MeV 

Harmonic  Oscillator  Length  Parameter  1.898  ~  .025  FM 


Single  Particle  Unit  Strength 


Slope 

-.4197  ~  . 

0173 

QEFF 

FF  ^ 

ERROR 

FF  CALE 

.594 

3.83E-3 

38.376 

3.01E-3 

.691 

3.66E-3 

24.494 

3.86E-3 

.792 

4.07E-3 

13.180 

4. 31E-3 

.894 

4. 35E-3 

15.276 

4.19E-3 

1.226 

1. 46E-3 

5.861 

1. 45E-3 

1.258 

1.18E-3 

5.693 

1.19E-3 

1.461 

2.02E-4 

23. 788 

1.72E-4 

1.072 

2.07E-5 

48. 809 

7.04E-6 

1.829 

8.27E-5 

19.814 

7.67E-5 

1.932 

1.13E-4 

14. 947 

l.lOE-4 

2.032 

1.20E-4 

13.470 

1.17E-4 

2.126 

9.87E-5 

12.452 

1.06E-4 

9.0081  -  .1023 

Chi  Square  1  Deg  of  Freedom  .3817 


Y 

P 

+  OR 

- 

.  3177 

.979 

.173 

.21 

.4300 

.798 

.0924 

.105 

.5648 

.741 

.0473 

.0506 

.7197 

.711 

.0524 

.0566 

1.3535 

.434 

.0125 

.0129 

1.4251 

.400 

.0112 

.0116 

1.9221 

.209 

.0236 

.0266 

2.5174 

.0967 

.0213 

.02  75 

3.0124 

-.274 

.026 

.0287 

3.3612 

-.417 

.031 

.0324 

3.7181 

-.568 

.037 

.0396 

4. 0701 

-.684 

.  0413 

.0440 

TABLE  V.2 


Level  Three,  =  2^  Excitation  Energy  3.499  MeV 

Harmonic  Oscillator  Length  Parameter  1.880  FM 
Single  Particle  Unit  Strength  1.811  “  .0783 


Slope 

-.3319  -  . 

0397 

Chi 

Square  1 

Deg  of 

Freedom  1 

.4224 

QEFF 

FF 

%  ERROR 

FF  CALC 

Y 

P 

+  OR 

.592 

1.53E-3 

74.835 

6.62E-4 

.  3097 

1.37 

.440 

T?80 

.694 

9.38E-4 

69.047 

8,90E-4 

.4256 

.  882 

.265 

.391 

.790 

1.09E-3 

8.210 

1.03E-3 

.  5515 

.  839 

.0338 

.0352 

.892 

9.59E-4 

7.401 

1.06E-3 

.  703 

.729 

.026 

.027 

1.221 

5.75E-4 

33.611 

5.33E-4 

1.317 

.  585 

.091 

.108 

1.253 

4.97E-4 

19.011 

4.67E-4 

1.  387 

.557 

.051 

.055 

1.456 

1.59E-4 

40.165 

1.47E-4 

1.873 

.  393 

.072 

.089 

1.667 

2.64E-5 

31.957 

1.73E-5 

2.4554 

.229 

.034 

.040 

1.825 

2.36E-6 

360. 343 

1.38E-7 

2.943 

.096 

.110 

- 

2.027 

1.48E-5 

33. 376 

3.77E-6 

3.640 

-.406 

.063 

.074 

2.122 

1.40E-5 

34.239 

5.28E-6 

3.978 

-.  523 

.083 

.098 

TABLE  V.3 


Level  Eight,  J 


Excitation  Energy  5.23  MeV 


QEFF 

FF 

%ERROR 

.588 

2.20E-4 

452.93 

.786 

1.62E-4 

194.17 

.889 

3.82E-4 

87.92 

1.247 

5.98E-5 

41.85 

1.45 

7.10E-5 

68.63 

1.818 

7.21E-6 

411.93 

1.921 

1.79E-5 

155.35 

2.115 

l.llE-4 

11.046 

TABLE  V.4 


Level 

Eleven,  J 

=  3 

Excitation 

Energy 

5.487  MeV 

Harmonic  oscillator  Length 

Parcuneter 

1.901 f 

.146fm 

Single 

Particle 

Unit  Strength  7.9467 

*  .511 

Slope 

-.2639  1  . 

1088 

Chi  Square/Deg  of 

Freedom 

.937 

QEFF 

FF 

%ERROR 

FF  CALC 

Y 

P 

+  OR 

.588 

8.25E-5 

1096.5 

2.63E-4 

.312 

.516 

1.27 

- 

.690 

7. 36E-4 

73.06 

4.96E-4 

.430 

.98 

.34 

.52 

.786 

8.72E-4 

35.51 

7.60E-4 

.558 

.91 

.  15 

.18 

.888 

7.28E-4 

45.31 

1.03E-3 

.712 

.68 

.  14 

.18 

1.214 

1.14E-3 

5.08 

1.13E-3  1 

.  3314 

.66 

.02 

.  02 

1.246 

l.lOE-3 

4.665 

1.07E-3  1 

.4025 

.64 

.02 

.  02 

1.  449 

5.41E-4 

9.017 

5.73E-4 

1.897 

.  48 

.02 

.02 

1.66 

1.59E-4 

10. 938 

1.72E-4 

2.  489 

.  33 

.01 

.02 

1.817 

6.11E-5 

30.661 

3.97E-5 

2.9805 

.26 

.04 

.  04 

1.92 

1.50E-5 

87.96 

8.52E-6 

3.3302 

.16 

.06 

.10 

2.02 

1.36E-5 

54.029 

2.72E-7 

3.6861 

.19 

.04 

.06 

TABLE  V.5 
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Level  Twelve,  J  =  2  Excitation  Energy  5.613  Kev 

Harmonic  Oscillator  Length  Parameter  2.21t  .178  FM 
Single  Particle  Unit  Strength  1.02272  .2438 


Slope 

-.411 

Chi  Square/Deg  of 

Freedom 

.  3187 

QEFF 

FF 

% ERROR 

FF  CALC 

Y 

P 

+  OR 

.786 

2.79E-4 

102.577 

3.66E-4 

.  7546 

.603 

.255 

- 

.888 

4.25E-4 

73.02 

2.95E-4 

-9631 

.727 

.229 

.349 

1.213 

2.14E-5 

128.25 

3. 34E-5 

1.797 

.209 

.107 

- 

1.245 

2.60E-5 

88.95 

2.19E-5 

1.8932 

.242 

.091 

.162 

1.659 

1.73E-5 

71.69 

9.46E-6 

3.362 

.515 

.160 

.24 

1.817 

3. 47E-6 

412.06 

9.99E-6 

4.032 

.  387 

.489 

- 

2.019 

4.08E-6 

144.789 

5.39E-6 

4.97 

.92 

.  514 

— 

TABLE  V.  6 
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Level  Thirteen,  J  =  4  Excitation  Energy  5.95  MeV 


Harmonic  Oscillator  Length 

Parameter 

1.737 

t  .126 

FM 

Single 

Particle 

Unit  Strength  5. 

9952  t  .3861 

Slope 

-.1766- 

.099 

Chi  Square/Deg  of 

Freedom 

.517 

QEFF 

FF 

%  ERROR 

FFCALC 

Y 

P 

+  OR 

_ 

.69 

1.67E-4 

649.39 

2.45E-5 

.359 

2.45 

4.26 

- 

1.212 

2.99E-4 

13. 862 

3.15E-4 

1.107 

.  784 

.05 

.05 

1.244 

3.44E-4 

10.466 

3.32E-4 

1.166 

.  809 

.04 

.04 

1.447 

3.90E-4 

12.03 

3.72E-4 

1.  578 

.739 

.04 

.04 

1.658 

2. 72E-4 

7.  076 

2.91E-4 

2.073 

.613 

.02 

.02 

1.816 

2.11E-4 

10.569 

1.91E-4 

2.487 

.588 

.03 

.03 

1.918 

1.29E-4 

12.378 

1.31E-4 

2.77e8 

.506 

.03 

.03 

2.018 

9.07E-5 

13.793 

8. 35E-5 

3.07 

.477 

.03 

.03 

2.113 

4. 45E-5 

15.144 

4.99E-5 

3.3665 

.383 

.028 

.03 
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are  lost  in  the  radiative  tails  and  in  the  background. 

Tables  V.T  through  V.9  show  the  experimental  form  factors  for  °Si 
extracted  in  the  present  analysis  and  the  form  factors  obtained  by  Whit- 
ner(5.^).  The  two  sets  compare  very  well.  Minor  differences  in  the  form 
factors  can  be  accomted  for  by  the  fact  that  Whitner  was  unable  to  fit  all 
levels  in  a  single  computer  run.  Hence,  at  higher  q  the  radiative  tails 
and  background  must  be  approximated.  The  comparison  indicates  that  no  ano¬ 
malies  have  crept  into  the  data  from  any  of  the  programs  used. 

It  is  known  the  ^^Si  is  not  a  perfectly  spherical  nucleus.  Therefore, 
the  ^^Si  states  are  most  probably  admixtures  of  various  spherical  shell 
model  components.  Comparing  the  extracted  polynomials  with  those  derived 
by  Donnelly  and  Haxton  places  experimental  limitations  on  the  possible  ad¬ 
mixtures  of  spherical  shell  model  states. 

In  general,  one  cannot  identify  uniquely  a  transition  from  the  shape 
of  its  fitted  polynomial.  In  particular,  any  first-order  polynomial  can  be 
generated  as  a  linear  combination  of  two  or  more  arbitrary  first-order  poly¬ 
nomials.  However,  if  the  fitted  polynomial  lies  close  to  one  of  the  derived 
polynomials  corresponding  to  a  particular  single-particle  transition,  there 
is  a  strong  suspicion  that  this  particular  single-particle  configuration 
contributes  significantly  to  the  transition.  An  accidental  fit  due  to 
particuleu-  combinations  of  other  single-particle  transitions  cannot  be  ruled 
out,  but  such  combinations  are  rather  unlikely.  It  is  in  this  spirit  that 
the  identifications  discussed  below  are  made. 

The  excited  state  in  ^^Si  is  j’*  *  2*.  The  fitted  polynomial  has  a 
.  + 

slope  of  -.42". 02  which  compares  to  a  slope  of  -.5  for  2S-ld.  The  lower  portion 
of  Figure  V.l  shows  the  fitted  polynomial  sis  well  as  2S-ld  and  Id-ld  transi- 


TABLE  V.7 


SILICON  28,0.0  MeV 


QEFF 

FF  EXP 

WHITNER  : 

.6079 

1.7741E-3 

- 

.7015 

2.0569E-1 

1.995E-1 

.7971 

1.179E-1 

1.147E-1 

.  8962 

6.0985E-2 

6.127E-2 

1.2346 

1.0492E-3 

1. 058E-3 

1.2668 

5.6132E-4 

5. 885E-4 

1.4698 

6.5732E-4 

- 

1.6804 

8.1412E-4 

8.222E-4 

1.8380 

5.024E-4 

4.598E-4 

1. 8945 

2.7193E-4 

- 

2.0402 

1.0750E-4 

1.119E-4 

2.1437 

3.7682E-5 

4.173E-5 

-fc-r 


TABLE  V.  8 

SILICON  28,1.78  MeV 


WHITNER  FF 


FF  EXP 

qeff 

2.260E-3 

- 

.604 

5.903E-3 

6.166E-3 

.698 

6.716E-3 

6.582E-3 

.794 

6.910E-3 

6.910E-3 

.896 

3.232E-3 

3.238E-3 

1.228 

2.714E-3 

2. 706E-3 

1.260 

7.217E-4 

6. 328E-4 

1.463 

3.424E-5 

3.624E-5 

1.674 

3.484E-5 

3.350E-5 

1.831 

7.947E-5 

7. 806E-5 

1.934 

1.044E-4 

1. 049E-4 

2.034 

1.077E-4 

1.090E-4 

2.128 

■» 


TABLE  V.  9 


SILICON  28,4.62  MeV 


QEFF  FF  EXP  WHITNER  FF 


.692 

4.33E-5 

4.100E-5 

.788 

6.716E-5 

7. 312E-5 

.890 

1.613E-4 

1.525E-4 

1.217 

4.164E-4 

4.256E-4 

1.250 

4.169E-4 

4. 468E-4 

1.453 

4.958E-4 

4.956E-4 

1.663 

3.590E-4 

3.593E-4 

1.821 

2.47E-4 

2.340E-4 

1.923 

1.650E-4 

1.528E-4 

2.023 

8.679E-5 

8.807E-5 

2.118 

4.665E-5 

4.984E-5 

tion  polynomials.  Also  Included  are  representative  polynomials  for  higher 
transitions,  It  is  important  to  note  that  the  form  factors,  at  least 

in  the  momentum  range  considered,  can  be  fitted  by  a  linear  polynomial. 

Hence,  the  higher  order  polynomials  of  the  transitions  do  not  contrib¬ 

ute.  This  particular  state  can  be  part  Id-ld  transition,  but  the  main  por¬ 
tion  is  most  likely  a  2S-ld  transition. 

I  *  ^ 

The  transition  to  the  second  excited  state  (J  =  2  )  appears  to  be  due 
more  to  Id-ld  transitions.  The  fitted  polynomial  has  a  slope  of  -.332  - 
.040  as  compared  to  a  slope  of  -.286  for  Id-ld  transitions.  Figure  V.2 
shows  tha  fitted  polynomial  as  well  as  the  derived  polynomials  for  2S-ld  and 
Id-ld  transitions. 

The  eleventh  state  has  j"  =  3  .  The  state  appears  to  be  a  Id-lf  trans¬ 
ition.  The  fitted  polynomial  is  a  straight  line  with  a  slope  of  -.264  - 
.109  as  compared  to  the  calculated  polynomial  slope  of  -.222.  Other  transi¬ 
tions  which  might  contribute  are  2S-lf,  2S-2p,  and  ld-2p  transitions,  but  the 
id-lf  transition  appears  to  be  dominant  (see  Figure  V.4).  Another  possibil¬ 
ity  for  creating  a  3~  state  is  a  Ip-ld  transition.  Such  a  possibility  is 
highly  unlikely  due  to  the  necessity  for  a  hole  buried  in  a  shell  below  the 
active  shell.  Furthermore,  this  transition  cannot  account  for  the  entire 
strength  because  its  polynomial  has  zero  slope,  and  as  indicated  above,  the 
fitted  polynomial  has  a  finite  negative  slope. 

The  twelth  state  is  ■  2^  and  appears  to  be  dominantly  a  2S-ld 

transition,  as  in  the  first  excited  state.  The  fitted  slope  is  -.411  - 
.104  whereas  the  slope  calculated  for  a  2S-ld  transition  is  -.5  (see  Figure 
V.5). 

x  -t- 

The  last  state  is  the  thirteenth  state  which  has  J  -  4  .  If  the 
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transition  were  purely  a  Id-ld  transition,  the  polynomial  would  be  a  hori¬ 
zontal  line.  However,  the  polynomial  has  a  slope  of  -.177  '  .099.  Since  the 
Id-ld  transition  does  not  contribute  to  this  slope,  only  higher  transitions 
can  be  responsible.  The  Id-lg  transition,  as  well  as  the  If-lf  transition, 
has  a  slope  of  -.182  and  either  or  both  could  possibly  be  responsible  (see 
Figure  V.6).  What  is  needed  are  data  taken  at  a  higher  momentum  transfer 
to  determine  where  the  minimum  occurs. 

Several  conclusions  can  be  reached  by  looking  at  the  2^  states,  the  3 
state,  and  the  4^  state,  separately.  The  2^  states  have  form  factors  which 
are  fitted  well  by  linear  polynomials.  Therefore,  the  single-particle  transi¬ 
tions  responsible  for  these  states  are  transitions  within  one  shell.  Transi¬ 
tions  of  2f>o  were  not  detected  within  this  momentum  range.  The  3  state 
appears  to  be  a  result  of  a  Id-lf  transition.  Again,  the  extracted  poly¬ 
nomial  is  linear.  The  4^  state  is  Interesting.  It  appears  as  if  one  of 
two  situations  is  occurring.  Either  there  is  some  If  strength  in  the  ^^Si 
ground  state  or  a  single-particle  transition  of  2KtJ  is  being  seen,  a  Id-lg 
transition.  Should  the  ground  state  have  some  If  strength,  then  a  If-lf 
transition  would  yield  a  polynomial  with  slope  -.182,  just  as  a  Id-lg  transi¬ 
tion  would.  Due  to  the  large  gap  between  the  Id  state  and  the  Ig  state,  it 
is  more  likely  that  the  ground  state  has  some  If  strength. 
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